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Using the latest quark Sivers functions extracted from the global analysis of available data on 
single transverse spin asymmetry (SSA), we calculate the SSA of Drell-Yan inclusive production of 
lepton pairs of invariant mass Q at both 4 < Q < 9 GeV and Q ~ Mz in p T p collisions at RHIC 
energies. We find that the features of the asymmetry An for Q ~ Mz are significantly different from 
that when 4 < Q < 9 GeV. The An near Z° pole is positive and sizable in the central rapidity region 
while the An at low Q is negative and only sizable in the forward rapidity region. We show that the 
size of An is sufficiently large for a wide range of Q around the Z° pole, and even with the lepton 
pair's invariant mass integrated from Q — 70 GeV to 110 GeV, the An is still close to 10 percent 
for the central and near forward rapidity region. We argue that the SSAs of Drell-Yan inclusive 
dilepton production at low Q and that near the Z° pole provide complementary information for 
testing the time-reversal modified universality (the sign change) of the Sivers functions. 

PACS numbers: 12.38.Bx, 12.39.St, 13.85. Qk, 13.88+e 



I. INTRODUCTION 

Single transverse-spin asymmetries (SSAs) in both semi-inclusive deeply inelastic scattering (SIDIS) and high 
energy hadronic collisions have attracted much attention from both experimental and theoretical sides in recent years. 
From the parity and time-reversal invariance of the strong interaction dynamics, the measured asymmetries in these 
collisions should be directly connected to the transverse motion of partons inside a polarized hadron. Understanding 
the dynamics behind the measured asymmetries should have the profound impact on our knowledge of Quantum 
Chromodynamics (QCD) and hadron structure 

Two QCD-based approaches for analyzing the observed novel SSAs, the higher twist coUinear factorization approach 
@> S EL 11 and the transverse momentum dependent (TMD) distribution approach [|| 0, H, H [Hj], have been 
proposed and bee n ap plied extensively in phenomenological studies. The first approach is based on the QCD coUinear 
factorization [l^, [Hi- By generalizing the successful leading power QCD coUinear factorization formalism to the 
next-to-leading power in the momentum transfer of the collision, non- vanishing SSAs could be generated by quantum 
interference between the real part and the imaginary part of the scattering amplitudes, and the SSA is proportional to 
twist-three transverse-spin-dependent multiparton correlation functions @, [H, 13 • This coUinear factorization based 
approach is more relevant to the SSAs of the processes with all observed momentum transfers Q 3> Aqcd- The 
other approach relies on the TMD factorization in QCD [15|, Ll6|, LUZI] - In this approach, the non- vanishing SSAs are 
attributed to the spin-dependent TMD parton distribution functions (PDFs), known as the Sivers functions [f|, or the 
TMD fragmentation functions, known as the Collins functions 0- The TMD factorization based approach is more 
suitable for evaluating the SSAs of the processes with two very different momentum transfers, Q\ Qi > Aqcd- 
The larger scale Q\ is necessary for using perturbative QCD while the lower scale Q2 makes the observable sensitive 
to the momentum of parton's transverse motion. These two approaches have their own kinematic domains of validity, 
which are complementary to each others. In the region where their validity overlaps, Q\ ^> Q2 Aqcd, these two 
seemingly different approaches should describe the same physics, and indeed, they were shown to be consistent for 
various processes where both factorizations were proved to be valid [18| . 

The predictive power of both approaches relies on the respective factorization. For the coUinear factorization based 
approach, the predictive power is a consequence of the infrared safety of all-order perturbatively calculatable partonic 
hard parts and the universality of all factorized and nonperturbative PDFs and correlation functions [12l . Il3j | . Unlike 
the collinearly factorized PDFs and correlation functions, the non-perturbative TMD PDFs or the Sivers functions 
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are not exactly universal, instead, they have the time-reversal modified universality [15j |. It was shown from the 
parity and time-reversal invariance of QCD that the Sivers functions in SIDIS and that in Drcll-Yan (DY) process 
should have the same functional form but an opposite sign - the time-reversal modified universality [la]. From the 
measurement of Sivers functions in SIDIS, the sign of Drcll-Yan SSA is uniquely predicted by the TMD factorization 
and the time-reversal modified universality of Sivers functions. 

The experimental test of the sign change of Sivers functions would provide a critical test for the TMD factorization 
approach and our understanding of the SSAs. Recently, HERMES and COMPASS experiments have performed 
measurements of Sivers functions in SIDIS [l9l |2C| . A new set of Sivers functions of various quark flavors was derived 
by the global analysis of available data [2l| . Future measurements of the SSAs in DY production have been planned 
(22I [23| . The SSAs of boson production in p^p collision at RHIC were also proposed to measure the DY-type 
Sivers functions (24[. Since it is difficult to reconstruct bosons by the current detectors at RHIC, we recently 
proposed to measure the SSA of single lepton decayed from W bosons at RHIC to test the time-reversal modified 
universality of Sivers functions [25| . We argued that with the large mass of W bosons M\y setting up the hard scale 
and the fact that the typical transverse momentum of W bosons |q_J produced at ^/s = 500 GeV at RHIC is much 
less than the W mass (|q_i_| <C M\y), the single lepton production from W boson decay at y/s = 500 GeV at RHIC is 
an ideal process for the TMD factorization. We demonstrated that the SSA of single lepton decayed from W bosons is 
sizable and its rapidity distribution is sensitive to the quark flavor of Sivers functions, and therefore, a good observable 
to test the sign change of the Sivers functions [25| . 

In this paper, we extend our previous work on the SSA of single lepton from the decay of W bosons to the SSA 
of inclusive lepton pair production with the pair's invariant mass near Z° boson mass Mz at *Js = 500 GeV at 
RHIC. In addition, we update the prediction of SSA of Drell-Yan lepton pair production with the pair's invariant 
mass Q between the J/ifj and T at RHIC energies by using the Sivers functions recently extracted from data of SIDIS 
experiments [2l| . We show that the asymmetry An changes sign from low Q to high Q in mid-rapidity region and is 
sizable. We find that the Am near Z° pole is significant for a wide range of Q around Mz, and the asymmetry could 
be used to test the time- reversal modified universality (the sign change) of the Sivers functions. 

The rest of the paper is organized as follows. In Sec. In] we derive the theoretical formalism for calculating the SSA 
of inclusive lepton pair production in hadronic collisions in the TMD factorization approach. In Sec. IIII1 we present 
our numerical predictions for the SSA of lepton pairs for RHIC kinematics. Finally, we summarize our findings, 
corresponding conclusions and the potential improvements for our calculations in Sec. IIVI 



II. SINGLE TRANSVERSE SPIN ASYMMETRY OF DILEPTON PRODUCTION 



The SSA of Drell-Yan lepton pair production via a virtual photon 7* has been studied previously in Refs. [22j, |23j, 
HEH2I- -"- n this paper, we study the SSA of inclusive lepton pair production via both a virtual photon 7*, a, Z° boson, 
and their interference, 



A^(p A , S) + B( PB ) -+W,Z -f] £+£-{Q) + X . 



(1) 



in hadronic collisions between a polarized proton A of momentum pa and transverse-spin vector S and an unpolarized 
proton B of momentum p B - The scale of hard collision Q in Eq. ([T]) is the invariant mass of the observed lepton pair. 
For testing the TMD factorization and the time-reversal modified universality of Sivers functions, we concentrate on 
the kinematic region where the transverse momentum of the produced lepton pair |qj_| is much less than the mass 
of the pair (|q^| <C Q), which is the region of phase space where the TMD factorization formalism is expected to be 
valid 

In terms of the TMD factorization, the leading order cross section for dilepton production can be written as 



dQ 2 dyd 2 qj_ 



1 



d 2 k a _L(i 2 k b± 5 2 (q ± - k a± - k bJ _)f q/A i (x a , k aJ _, S)f q / B (x h , k b ±) <r, 
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(2) 



where runs over all light (anti)quark flavors, y is the rapidity of the observed lepton pair, and kb± = |kf,jj is 
the magnitude of parton transverse momentum. The lowest order partonic cross section oyj 9 in Eq. @ receives the 
contribution from both the virtual photon and the Z° boson, as well as their interference, and is given by, 



1 

n7s 



(3) 



where s = (pa +Pb) 2 , ol and Gf are the electro-magnetic and Fermi weak coupling constants, respectively, e q is the 
quark fractional charge of flavor q, N c = 3 is the number of color, v q ^ and a q (i) are the vector and axial couplings of 
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the Z° boson to the quark (lepton), respectively, and 

R 



Ml 



Q 2 -M 2 +iM z T z 



(4) 



with Mz and Y z the mass and decay width of the Z° boson, respectively. The parton momentum fractions in Eq. ([2]) 
are given by, 



Q 



Xb 



(5) 



to the leading power in q\/Q 2 . 

Following the notation of Ref . 21| , we expand the spin dependent TMD quark distribution as, 

fq/Al( x >k±,S) = f q /A(x,k±) + -A N f q/AX (x,k±)S ■ (p A x k±), 



(6) 



where = |kj_|, p A and kj^ are the unit vectors of pa and k^, respectively, fq/A{ x ik±) is the spin-averaged TMD 
quark distribution function of flavor q, and A N f q /Ar(x, k±) is the Sivers function of flavor q. Substituting Eq. ^ to 
Eq. |2]), we obtain for the spin-dependent cross section as 



dAa A 



(S) 



dQ 2 dy d 2 qj 



da A iB^£+e-+x(S) da A 



t b-*i+i-+x 



(-S) 



dQ 2 dy d 2 q± dQ 2 dy d 2 q± 

J d 2 k a± d 2 k b± S 2 (q ± - k a± - k bJ _) S • (p A x k Q 



'xA N f q/A r (x a ,k a ±) fq/ B {xb, h±) a% q , 
and corresponding spin-averaged cross section as 

dcf AB^l+t~+X _ 
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= ^2 J d 2 k a xd 2 k b ±5 2 (qx -k aJ _ - k b ^) f q / A {x a ,k al _) f q/B (x b ,k bl _)&l q 



dQ 2 dy d 2 q± 

From Eqs. ([7]) and ©, we obtain the conventionally defined SSA as 

dAa A T B^i+e-+x(S) / da A B^e+e-+x 



A 



N 



dQ 2 dy d 2 q\ 

The commonly used weighted asymmetry is defined as 



dQ 2 dy <i 2 qj 



A 



sin(</> — </> s 
N 



dAcr A T B ^e+i- + x(S) 
dQ 2 dy d 2 q± 



da 



AB^i+£-+X 



dQ 2 dyd 2 q_ 



± 



(7) 



(8) 



(9) 



(10) 



where (f> is the azimuthal angle of the lepton pair of transverse momentum qj_ and <j) a is the azimuthal angle of the 
spin vector S of incoming polarized hadron. 

In order to evaluate the SSAs in Eqs. © and (jTUJ) for inclusive lepton pair production in p^p collisions, we adopt 
the parameterization of TMD parton distributions introduced in Ref. [2l[ , 



f q/h {x,kx) = f q (x) 



1 



7r(fc] 



3 -fci/<fci> 



> N fg/hf S &k±) = 2M q {x)h{k 1 _)f q/h {x,k 1 _l 

h(k x ) = v^i— e~ fc °/ Ml 
V ' Mi 



(11) 
(12) 
(13) 



where f q {x) is the standard unpolarized parton distribution of flavor q, (kj_) and Mi are fitting parameters, and Af q (x) 
is a fitted distribution given in Ref. [2lj . 

With the Gaussian ansatz for the transverse momentum dependence of the TMD distributions, one could carry out 
the integration d 2 k a ±d 2 k b ± in Eqs. © and © analytically, and obtain 



A N = S • (p A x qj_)- 



f2e 2(k 2 



M x [(k\) + (k 2 )} 2 



p-jEg^f [-Af q (x a )]f q / A (x a )fq/ B (x b ) 



J2q 6l q fq/A(x a )f q/B (x b ) 



(14) 
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where (k^) = M\ (fcj_)/[Mf + ancl the "— " sign in front of N q (x a ) is due to the sign difference between the 

Sivers function in DY process and that in SIDIS process. 

For our numerical results in next section, we choose a frame in which the polarized hadron A moves in the +z- 
direction, and S and qj^ along y- and x-directions, respectively. It is important to realize that the weighted asymmetry 
defined in Eq. (fl~0)) differ from An in Eq. ([9]) by a minus sign, 

A p4>-4>s) = _ An (15) 

in our choice of frame. 



III. PHENOMENOLOGY 

In this section, we present our numerical estimates for the SSA of inclusive lepton pair production in p^p collisions 
at RHIC energies by using Eq. (fl4|) . We use GRV98LO parton distribution [28[ to be consistent with the usage of the 
TMD distributions of Ref. [2l|. 

In Fig. [I] we plot the asymmetry An for the lepton pair production with the pair's transverse momentum integrated 
up to qt = |q_i_| = 1 GeV at RHIC energy y/s = 200 GeV (solid lines) and — 500 GeV (dashed lines). In the figure 
on the left, we plot An as a function of rapidity y of the lepton pair with the pair's invariant mass Q integrated from 
4 GeV to 9 GeV. Similarly, in the figure on the right, we plot An as a function of the pair's invariant Q with the 
pair's rapidity integrated from to 3. Since the mass of Z° boson is so much larger than the lepton pair's invariant 
mass Q considered in Fig. [TJ the An at yfs = 200 GeV (solid lines in two figures) are effectively the same as those 
presented in Ref. (23[, in which only the virtual photon channel was considered. That is, as expected, the Z° boson's 
contribution to the production rate as well as the asymmetry is not important for the region where Q <JC Mz- 




y q 

FIG. 1: SSA of lepton pair production as a function of the pair's rapidity y (left) and invariant mass Q (right) at y/s = 200 
GeV (solid lines) and y^s = 500 GeV (dashed lines). We have integrated the pair's transverse momentum from to 1 GeV. 

As shown in Fig. [I] An of inclusive lepton pair production in this small Q region is negative. Unlike the TV ± -boson 
production (or the lepton production from the decay of VF-boson) (25| . Drell-Yan inclusive lepton pair production 
via a virtual photon or a Z° boson is not very good in separating contributions from different quark or antiquark 
flavors. The SSA in this low Q region is dominated by the virtual photon channel and is proportional to a sum of 
quark (antiquark) Sivers function of all flavors weighted by the quark fractional charge square and corresponding 
spin-averaged antiquark (quark) distribution. Since the it-quark and <i-quark Sivers functions have the opposite sign 
[2l| , the sign and size of the SSA of inclusive dilepton production is a result of an incomplete cancelation between 
contributions from different quark flavors. For the latest parameterization of the Sivers functions obtained in Ref. [21|, 
the it-quark Sivers function for Drell-Yan type processes is negative while d-quark Sivers function is positive, and the 
absolute size of it-quark Sivers function is slightly smaller than <i-quark Sivers function. Since the contribution to the 
SSA from it-quark Sivers function is weighted by a larger fractional charge square, it wins over the contribution from 
ci-quark Sivers function, and consequently, we have a negative An as shown in Fig. [TJ The rapidity dependence of 
An is a consequence of the x-dependence of the Sivers functions extracted from the global fitting in Ref. [2l|. The 
measurement of the rapidity dependence of the SSA of Drell-Yan lepton pair production could provide information 
on both the sign and the functional form of the Sivers functions to test the time-reversal modified universality of the 
Sivers functions. 
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In Fig. [T] we also plot the SSA of Drell-Yan lepton pair production at y/s = 500 GeV (dashed lines). It is clear 
that the asymmetries at y/s = 500 GeV carry the same features as those at y/s — 200 GeV. However, due to the 
larger collision energy while keeping the same range of qx and Q, the same rapidity y corresponds to different values 
of momentum fractions of colliding partons. That is, the rapidity dependence of the SSA at different collision energy 
provides additional information on the functional form of the Sivers functions. 

The SSA of Drell-Yan lepton pair production is small in the central region of the collision (y ~ 0) due to the 
cancelation of the contributions from different quark flavors. The asymmetry becomes sizable for the experimental 
measurement only in the forward region (y ~ 1 — 3), as shown in Fig. [U where the production rate for the lepton pair 
is significantly reduced. Since the hard scale Q is not very large in this low Q region, the theoretical uncertainty for 
(or correction to) the leading power TMD factorization formalism becomes larger in the large y forward region. In 
principle, the TMD factorization works better when the hard scale Q is much larger than qx- But, in practice, the rate 
for lepton pair production quickly diminishes when the invariant mass Q and/or rapidity y increase while the collision 
energy y/s is fixed. However, with the new collision energy at y/s — 500 GeV, RHIC can produce sufficient numbers 
of W ± and Z° bosons. The large mass of and Z° makes the corrections to the TMD factorization formalism 
smaller and the theoretical calculations more reliable. The SSA of inclusive lepton pair production at Q = Mz, which 
is mainly from the decay of Z° , could be a very good complementary observable to the SSA of IV ± production to test 
the TMD factorization in QCD and the time-reversal modified universality of the Sivers functions. 

In Fig. [21 we plot the SSA of lepton pair production at Z° pole (Q — Mz) at RHIC at y/s = 500 GeV as a function 
of the pair's transverse momentum q? (left) and the rapidity y (right). We find that the SSA of inclusive lepton pair 
production at Z° pole is positive, in contrast to the negative SSA in Fig. [T] Different from the low Q region, the 
production of lepton pair at Q = Mz is dominated by the Z° channel. Although both u- and <i-quark Sivers function 
in the Z° channel contribute to the SSA in the same way as that in the virtual photon channel, the relative weight 
between their contributions is different. Instead of the fractional charge square e 2 in the virtual photon channel, the 
contribution of quark flavor q in the Z° channel is weighted by v 2 + a 2 with 

1 4 . 2 1 

v u = w a " = 2' 

Vd = -^ + ^sin 2 6 W a d = -i. (16) 

Unlike the virtual photon channel, where > e 2 ,, the relative weight for the contribution of it- and <i-quark Sivers 
function in the Z° channel is opposite, v^+a^ < v\ + a 2 ,. Since the absolute value of d-quark Sivers function is larger 
than that of it-quark Sivers function as found in Ref. [2l| , the contribution from the d-quark Sivers function wins over 
that from it-quark Sivers function in the Z° channel, and consequently, An > when the lepton pair's invariant mass 
is near the Z° pole, Q ~ Mz, as shown in Fig. O 
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FIG. 2: SSA of lepton pair production at Z° pole Q = Mz as a function of the pair's transverse momentum qr (left) and 
rapidity y at y/s = 500 GeV. 

As shown in Fig. [2 An of inclusive lepton pair production with the pair's invariant mass Q = Mz in p^p collisions 
at y/s = 500 GeV is very significant at qx ~ 1 GeV in the central y region. The absolute size of the asymmetry in 
the central region is much larger than those shown in Fig. [1] for the low Q region. The SSA of inclusive lepton pair 
production near Z° pole could be an equally good observable as the low Q Drell-Yan lepton pair production to test 
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the time-reversal modified universality of the Sivers functions, if there is a sufficient rate of producing Z° bosons at 
RHIC. 

To better understand the sign change of An of inclusive lepton pair production from the low Q region to the region 
near Z° pole, we plot the An as a function of Q in Fig. [3J (left). The figure clearly demonstrates the transition from 
the negative An at a low invariant mass Q to the positive An when Q is near Z° pole. In our calculation, the lowest 
order partonic cross section, <5q 9 in Eq. ([3]), determines the relative weight between the contributions of it-quark and 
d-quark Sivers function at a given dilepton invariant mass Q. 




q y 

FIG. 3: Left: SSA of lepton pair production as a function of the pair's invariant mass Q. Right: SSA of lepton pair accumulated 
around Z° pole as a function of rapidity y. 

From Fig. [3J (left), the ^4^ is sufficiently large even when the pair's invariant mass is not exactly sitting at the Z° 
pole. Therefore, we could have a relatively wide bin of Q to improve the event rate of inclusive lepton pair production 
at large Q. In Fig. [3J (right), we plot the SSA as a function of rapidity y for all lepton pairs whose invariant mass Q 
are between 70 GeV to 110 GeV. It is clear that the SSA of inclusive lepton pair production at y/s — 500 GeV for a 
wide bin of invariant mass Q around the Z° pole is close to 10% at the central rapidity, and could be measurable in 
the future at RHIC. 



IV. CONCLUSION 

We studied, in terms of TMD factorization formalism, the single transverse spin asymmetry of inclusive lepton pair 
production in p^p collisions at RHIC energies for the kinematic region where the lepton pair's transverse momentum is 
much less than the pair's invariant mass, qx *C Q, the condition that is required for justifying the TMD factorization. 
In our calculation of the SSA, we included contributions to the lepton pair production from both the virtual photon 
and Z° boson channel as well as the contribution from their interference. 

By using the latest Sivers functions extracted from the global analysis of available data [2l| , we evaluated the SSA 
of Drell-Yan lepton pair production for 4 < Q < 9 GeV at y/s = 200 GeV and 500 GeV. The asymmetry An at 
this low Q region is close to zero for the central and negative rapidity region due to the cancelation between the 
contribution from it-quark and d-quark Sivers functions. The asymmetry An is negative and sizable in the forward 
rapidity region (y ~ 1 — 3) [23l |. 

We also calculated the SSA for inclusive production of lepton pairs with the pairs' invariant mass near Z° boson 
pole in p^p collisions at y/s — 500 GeV. We found that the features of the asymmetry An for Q ~ Mz are significantly 
different from that when 4 < Q < 9 GeV. As shown in Figs. Q] and [2 the asymmetry An near Z° pole is positive 
and sizable in the central rapidity region while the A n at low Q is negative and only sizable in the forward rapidity 
region. The sign change of the asymmetry is due to the transition of contributions from the virtual photon channel 
dominance to Z° channel dominance when Q increases, which leads to the change of the relative weight of the 
contributions from it-quark and (i-quark Sivers function that have opposite sign. Since the effective weight of the 
different flavor contributions to the asymmetry varies when the hard scale Q changes, we conclude that the SSAs 
at low Q region and that near Z° pole provide complementary information on both the sign and the functional form 
of the Sivers functions as well as their flavor dependence. 

We noticed that the size of asymmetry An for inclusive lepton pair production is sufficiently large for a wide range 
of Q around the Z° pole. We found that the An with lepton pair's invariant mass integrated from Q = 70 GeV to 
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110 GeV is close to 10 percent for the central and near forward rapidity region. This large bin width for the lepton 
pair's invariant mass helps to increase the event rate, and the size of should be measurable in the future at RHIC. 
The successful measurement of SSA of inclusive lepton pair production could provide a crucial test of the time-reversal 
modified universality (or the sign change) of Sivers functions and the validity of TMD factorization approach to the 
novel phenomenon of SSA. 

Our results and conclusions for the SSA of inclusive lepton pair production rely on the TMD factorization and 
the accuracy of the extracted spin-averaged TMD parton distributions and the Sivers functions as well as their 
factorization scale dependence. Unlike the parton distribution functions in the collinear factorization approach, whose 
factorization scale dependence is well-understood and determined by the DGLAP evolution equation, the factorization 
scale de pen dence of TMD parton distributions is not yet fully understood. Without reliable theoretical input, authors 
of Ref. [2l[ assumed that the Sivers functions obey the spin-averaged DGLAP evolution equation in their global 
analysis of extracting TMD parton distributions. Although it is difficult to judge how good the assumption is, we 
noticed that the evolution equation for the twist-3 quark-gluon correlation functions relevant to the SSA in collinear 
factorization approach, which are one-to-one related to the first moment of the quark Sivers functions , is different 
from the spin-averaged DGLAP evolution equation [25|,[3(|. But, the difference is relatively small as long as the parton 
momentum fraction x is not too small [25j . Therefore, we think that using the Sivers functions from Ref. (2lj , which 
fit all existing data, is a reasonable choice that we can make now, and the general features of the SSA of inclusive 
lepton pair production presented in this paper should be valid. For the precision test of TMD factorization in QCD, 
it is of course very important to understand the factorization scale dependence of the Sivers functions or TMD parton 
distributions in general, which is beyond the scope of this paper. 

When the hard scale of the collision, Q, the invariant mass of the produced lepton pair is so much larger than the 
observed transverse momentum of the pair, Q ^> qt, the Q and qx dependence of the production rate is dominated 
by the gluon shower right before the hard quark-antiquark annihilation to produce the lepton pair. The Q and qr 
dependence of the lepton pair production in this kinematic region is dominated by the hi 2 (Q 2 / q^)-ty~pe of leading 
logarithmic contribution and can be calculated by using the Collins-Soper-Sterman (CSS) resummation formalism in 
the 6-space, the Fourier transform of the qr space [31| . Although the Fourier transformation from 6-space calculation 
to the qr-space could be sensitive to the non-perturbative physics at large b, the resummation formalism does have an 
excellent predictive power if the Fourier transform of the 6-space distribution is dominated by the small b region where 
the perturbative QCD calculation is effectively done in collinear factorization and reliable [32J. The CSS resummation 
formalism should work well for processes with a large hard scale Q and large enough phase space for the gluon shower 
[33| . The formalism has been successful in interpreting the data on and Z° production in spin-averaged hadronic 
collisions jolj . The generalization of the resummation formalism to the spin-dependent cross section with transverse- 
spin is nontrivial when the partonic hard part is not azimuthal symmetric, but, it should be a very interesting and 
important project for the future publication [35l . |36| . 
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